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Abstract 

We propose a new framework for the study of continuous time dynamical systems on net- 
works. We view such dynamical systems as collections of interacting control systems. We show 
that a class of maps between graphs called graph fibrations give rise to maps between dynamical 
systems on networks. This allows us to produce conjugacy between dynamical systems out of 
combinatorial data. In particular we show that surjective graph fibrations lead to synchrony 
subspaces in networks. The injective graph fibrations, on the other hand, give rise to surjective 
maps from large dynamical systems to smaller ones. One can view these surjections as a kind 
of "fast / slow" variable decompositions or as "abstractions" in the computer science sense of the 
word. 
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1 Introduction 
1.1 Overview 

A fundamental question in the study of dynamical systems is to determine the existence and 
properties of a map that intertwines the dynamics of two different systems. Stated concretely, 
given two manifolds M, N, and two flows ipt : M — )• M, ipt'-^ does there exist a map 

h: M N such that 

h o ift = ij^t ° hi {1.1.1) 

Equivalently, given two manifolds and two vector fields X: M — )■ TA/, Y : N ^ TN, does there 
exist a map h such that 

dhoX = Yohl (1.1.2) 

If there is an h that satisfies (jl.l.ip or (jl.l.2p , then we call this a map of dynamical systems. Given 
such a map /i, we would like to understand its properties and to compute it explicitly. 

A common restriction requires that h be invertible. In this case, it is said that we have exhibited 
a conjugacy between the two dynamical systems, and this means that all of the dynamical features 
of the flow are the sam^. The notion of conjugacy of dynamical systems goes back at least to 
Poincare [IHS], it was further developed by Smale and collaborators (HE], and is now the basic 
notion in modern dynamical systems theory. 

A more general notion of conjugacy arises from the relaxation of the assumption of invertibility; 
here, the existence of the map h still produces significant information. For example, the flow ipt 
has a fixed point iff we can exhibit a map h: {*} — > A^, where {*} is a one-point set, and h 
satisfies (11.1. ID . If M = and the flow tpt is given by ipt{e^^) = e^Trit/Tgie^ then the existence of 
h: M N satisfying (jl.l.ip amounts to the flow ipt having a periodic orbit of period T. 

It is also common for h to be chosen to be surjective. In this case the map is typically termed 
a semiconjugacy [2l[8] and certain nice properties follow. We do not expand on this here, but we 
will exploit the existence of semi-conjugacies for networked systems below in Section O 

The question of determining whether a map relating two dynamical systems exists, and what 
its properties might be, is exceedingly difficult [9l412j in general. In some cases, even if such an h 
is known to exist, determining its form (or even qualitative properties) can be challenging. 

In a different direction, dynamical systems defined on networks have become the fundamental 
object of study across a variety of fields. Some examples include the design of communications 
networks [13]; cognitive science, computational neuroscience, and robotics (see, for example |14f - 

^Note that what we mean by "the same" depends on the category in which we work. For instance if /i is a 
homeomorphism then we say that the dynamical systems are topologicaUy the same. Many of the imphcations of the 
existence of a conjugacy are worked out in [T}{3|- 
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I19j): gene regulatory networks |20H22j and more general complex biochemical networks [23]; and 
finally in complex active media |24fl28j . 

There are multiple definitions in the literature of what it means to define a "dynamical system on 
a network" and we will not compare them here. A common thread running through these definitions 
is that imposing a network structure on a dynamical system should mean that component j of the 
system depends upon component i of the system iff the underlying graph has an edge i ^ j- 

In this paper we show that an imposition of a "network structure" on a dynamical system allows 
us to produce maps between dynamical systems in a precise, computable and combinatorial manner 
from finite data. Thus the purpose of this manuscript is twofold: first, to present a notion of a 
dynamical system "consistent with a graph"; second, to show that certain maps between graphs 
induce maps between the dynamical systems that live on them. In particular, we will show below 
that all graph maps that respect a particular combinatorial structure induce maps between the 
dynamical systems living on these graphs. 

We focus on the case where the dynamics is modeled by vector fields on manifolds. The 
interactions of subsystems are coded by directed (multi-) graphs with "labels". These labels in 
particular, assign to each node of a graph the phase space of the relevant subsystem. The ideas of 
the paper can be extended to both discrete-time, hybrid and stochastic systems, and we plan to do 
so in future work. 

As stated above, the main result of the paper is the construction of maps of dynamical systems 
from maps of labeled graphs. In particular we show that surjective maps of graphs, such as the one 
arising from quotienting a graph by an appropriate equivalence relation, give rise to embeddings 
of dynamical systems; second, injective maps of graphs give rise submersions of the corresponding 
phase spaces and surjective maps of dynamical systems. The former is very useful in character- 
izing the "modularity" of a networked dynamical system; the latter gives a precise mathematical 
formulation of some intuitive notions of whether and how we can think of a large dynamical system 
driven by a subsystem. 

1.2 Background and previous work 

The present paper is inspired by several distinct bodies of work that are well known in the applied 
mathematics communities. 

The first body of work has been mainly applied to chemical reaction systems, and, in some 
specific cases, to Petri nets; this work has been used in both the deterministic and stochastic 
settings. One of the earliest results in this direction is the "zero deficiency theorem" of Feinberg [291 - 
[34], first used to show the existence of stable equilibria in biochemical systems and then expanded to 
statements about the existence and structure of the equilibria in stochastic biochemical systems [35] . 
This type of methodology has also been expanded to Petri nets |36fl38j and models describing gene 
regulatory networks [21| . 

The second body of work is due to Golubitsky, Stewart, and various collaborators |39H64j . These 
authors considered a notion of ODEs (vector fields defined on Euclidean spaces) that were consistent 
with a graph structure. The resulting networks are called coupled cell systems and the approach the 
groupoid formalism. The main idea was to consider "balanced" equivalence relations on the vertices 
of a graph. They showed that these relations lead to the existence of certain invariant subspaces 
termed "poly diagonals." We will see that the quotient maps resulting from balanced equivalence 
relations are instances of graph fibration in the sense of Boldi and Vigna [65J. However, while we 
are greatly indebted to this body of work for its intellectual inspiration, we also point out that 
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our approach differs from this work in some very specific ways. We elucidate the connections and 
contrasts in Remark 14.3.31 below. 

An alternative approach to coupled cell systems has been developed by Field and collabora- 
tors |66l - l68j . They also considered ODEs and other types of dynamical systems consistent with 
directed graphs. This approach is considered broadly equivalent to that of Golubitsky et al. 

1.3 The contributions of the paper 

In this paper we propose a new framework for the study of continuous time dynamical systems on 
networks. We view such dynamical systems as collections of interacting control systems. We show 
that a class of maps between graphs called graph fibrations give rise to maps between dynamical 
systems on networks. This allows us to produce conjugacy between dynamical systems out of 
combinatorial data. While the current work is certainly inspired by the methods and results of 
both the Feinberg et al., Golubitsky et al. and, to a lesser extent, Field et al. groups, our approach 
and results differ in several important respects. 

1. Our basic philosophy is that of category theory — so rather than study dynamical systems 
one at a time we aim to study maps between all relevant dynamical systems at once. To 
quote Silverman [69]: 

A meta-mathematical principle is that one first studies (isomorphism classes 
of) objects, then one studies the maps between objects that preserve the objects' 
properties, then the maps themselves become objects for study and one tries to 
put a "nice" structure on the collection of maps (often modulo some equivalence 
relation) . 

2. Our set-up is coordinate-free and works for vector fields on manifolds, not just M". In this 
case we are enlarging both the scope of the Feinberg et al. work (polynomial vector fields 
on positive orthants) and that of Golubitsky et al. (vector fields on Euclidean spaces). This 
aspect of our approach is similar to the work of Field et al. There are several motivations for 
working on manifolds as opposed to ODEs living in Euclidean spaces. These include dealing 
effectively with constraints, and extending the results to the setting of geometrical mechanics. 

3. It will be evident from the construction below that the quotient maps of graphs by balanced 
equivalence relations of |61j are special cases of graph fibrations — they are the surjective 
graph fibrations. However, even in the case of surjective graph fibrations our maps of dynam- 
ical systems have the opposite direction from the maps in the groupoid formalism. Rather 
than restricting from polydiagonals we extend from polydiagonals. This allows us to deal 
with surjective and general graph fibration on the same footing. 

1.4 Motivating example 

Consider an ODE in (R")^ of the form 

xi = f(x2), X2 = f{xi), x-i = f{x2) (1-4.1) 
for some smooth function / : M" M". That is, consider the fiow of the vector field 
F : (M")3 ^ (M")3, F{xi,X2,xs) = (/(xa), /(xi), /(xs)). 
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It is easy to check that F is tangent to the diagonal 



M" ~ A = {(xi,rE2,rE3) G {W^f \ xi = = X3} 
and that the restriction of the flow of F to A is the flow of the ODE 

u = f{u). 

One can also see another invariant submanifold of F: 



~ A' = {{xi,X2,x-i) G (M")3 I 
On A' the flow of F is the flow of the ODE 

Moreover the projection 



Xi = X3}. 



TT 



pn\3 



f A', 1t{xi,X2,X3) = {xi,X2, Xi) 



intertwines the flows of F on and on A'. We have thus observed two subsystems of 

and three maps between the three dynamical systems: 



n\3 



(A,F| 



((M")3,F) 



(A',F|a' 



(1.4.2) 

Where do these subsystems and maps come from? There is no obvious symmetry of (M")^ that 
preserves the vector field F and fixes the diagonal A and thus could account for the existence of 
this invariant submanifold. Nor is there any F-preserving symmetry that fixes A'. In fact the 
vector field F does not seem to have any symmetry. The graph G recording the interdependence 
of the variables {xi,X2,X3) in the ODE (jl.4.ip has three vertices and three arrows: 



G 




(1.4.3) 



The graph has no non-trivial symmetries. Nonetheless, the existence of the subsystems (A,F|a), 
(A',F|a') and the whole diagram of the dynamical systems (jl.4.2p can be deduced from certain 
properties of the graph G. There are two surjective maps of graphs: 



:G^aO 



and 



ip : G 




with ^ defined on the vertices by ip{2) = b, V'(l) = a = and one embedding 





We can collect all of these maps into one diagram 

A comparison of (jl.4.2p and (jl.4.4p evokes a pattern: for every map which intertwines dynam- 
ical systems in (ll.4.2p . there is a corresponding map of graphs in (I1.4.4P with the arrows reversed, 
and vice versa. 

The same pattern holds when we replace the vector space by an arbitrary manifold M. 
Given a pair of manifolds U and A'^, we think of a map X : U x N ^ TN with X{u, n) € TnN as 
a control system with the points of U controlling the the dynamics on N. Now consider a vector 
field 

F ^ T{M^) = TM x TM x TM 

of the form 

F(xi,X2,X3) = (/(X2,3;i),/(X1,X2),/(X2,X3)) 

for some control system 

f : M X M ^ TM, with f{u, v) e T^M. 

Then once again the three maps of graphs in the diagram (I1.4.4P give rise to maps of dynamical 
systems 

{Am,F\^Jc ^{M\F) ,' {A\„F\^,) (1.4.5) 

What accounts for the patterns we have seen? Notice that the dynamical systems ()1.4.5p are 
constructed out of one control system f : M x M TM. At the same time, in each of the graphs 
in ()1.4.4p . every vertex has exactly one incoming arc. This is not a coincidence. The rough idea for 
the technology which generalizes this example is this: if we have a dynamical system made up of 
repeated control system modules whose couplings are encoded in graphs, then the appropriate maps 
of graphs lift to maps of dynamical systems. Making this precise requires a number of constructions 
and theorems; these make up the bulk of this paper. 

1.5 Main ideas of the paper 

We study dependency and modularity of networks and their effect on the concomitant dynamical 
systems. 

By dependency we mean the following. Each node in a network corresponds to a single dynamical 
variable living on a particular manifold. We then require that the variable corresponding to node 
i can depend on the variable corresponding to node j iff there is an edge in the graph from node 
j to node i. We give a more precise description of this requirement below and we will denote the 
space all vector fields with this property by S{G,V); see Section [2.41 below for a definition. 

The rough idea of modularity is that if we ever have multiple nodes of the graph that are "the 
same" and have "the same" inputs, then we require that these nodes are interchangeable in the 
dynamical system. Speaking more precisely, we will assume that in each network, each node has a 
"type" (it will, in fact, be a manifold attached to this node which corresponds to the phase space 
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of the variables associated to that node), and if we ever see two nodes, each with type x, with 
n inputs in the graph, such that these inputs are of type xi, . . . ,Xn, then the vector field defined 
on these two nodes must depend on their inputs in exactly the same manner. We will denote all 
vector fields that respect this principle of modularity as V(G, V) and give a precise definition of 
these vector fields in Section [3l 

In the example in the previous subsection, the principle of dependency tells us that the system 
living on the graph in (jl.4.3p must be of the form 

Xl = fixi,X2), X2=g{x2,Xi), X3 = h{x3,X2), (1.5.1) 

but xi cannot depend on X3, for example. The principle of modularity tells us that the functions 
/, g, h must all be the same, i.e. that 

Xl = f{xi,X2), X2=f{x2,Xi), X3 = f{x3,X2), (1.5.2) 

Therefore, all systems in S(G,P) must satisfy (|1.5.ip . but those in Y{G,'P) must satisfy the stricter 
requirements (ll.5.2p . 

2 Networks and dynamics on networks 

The goal of this section is to define networks and dynamics on networks in the context of continuous 
time dynamical systems. It is not uncommon to read that "a network is a graph." This could not 
be a complete story since by a network one usually means a collection of interconnected subsystems. 
We now make this precise. 

2.1 Graphs, manifolds and networks 

Throughout the paper graphs are directed multigraphs, possibly with loops and multiple edges 
between nodes. More precisely, we use the following definition: 

2.1.1 Definition. A graph G consists of two sets Gi (of arrows, or edges). Go (of nodes, or vertices) 
and two maps s, i : Gi — )■ Go (source, target): 

G = {s,t:Gi^Go}. 

We write G = {Gi ^ Go}. A graph G is finite if it has finitely many arrows and edges. 

2.1.2 Definition. A map of graphs ip : A B from a graph A to a graph i? is a pair of maps 
ifi : Ai ^ Bi, ipo : Aq ^ Bq taking edges of A to edges of B, nodes of A to nodes of B so that for 
any edge 7 of A we have 

^o{s{7)) = s{ipi{j)) and (po{t{j)) = t{ipi{j)). 
We will usually omit the indices and 1 and write 93(7) for fii'j) and ^{a) for ipoia). 

2.1.3 Remark. The collection of graphs and maps of graphs form a category Graph. The subcol- 
lection of finite graphs and maps of graphs forms a full subcategory FinGraph. 
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To construct a network from a graph we need to attach phase spaces to its vertices. Since we are 
interested in continuous time dynamical systems, we choose phase spaces to be (finite dimensional 
paracompact Hausdorff) manifolds. Other choices, of course, may also be reasonable, such as 
coordinate vector spaces or manifolds with corners. 

2.1.4 Definition (Network). A network is a pair (G, V) where G is a finite graph and "P is a function 
that assigns to each node a G Go of G a manifold V{a). We refer to 7^ as a phase space function. 
Note that the target of the function V is the collection Man of all manifolds: V : Go — > Man. 

A map of networks from (G, V) to (G', V') is a map of graphs : G — )• G' so that 

V' oip = V. 

2.1.5 Remark. It is easy to see that composition of two maps of networks is again a map of 
networks. In other words networks form a category. We denote it by FinGraph/Man. 

Given a network (G, V) as defined above, a state of the network is completely determined by 
the states of its nodes. Hence the total phase space of the network should be the product 

P(G,P) := fl V{a). 

aeGo 

Note, however, a small issue: there is no natural ordering on the vertices of the graph G. We could 
choose an ordering (ai, . . . ,0^) of the vertices and define the total phase space as the Cartesian 
product 

n 

P(G,P) :=llV{a,). 
1=1 

However, it will be convenient not to choose an ordering of vertices and use a slightly different 
notion of product. This version of the product is used, for example, in chemical reaction network 
literature |37| . 

2.1.6 Definition. Given a family {Ms}s£S of manifolds indexed by a finite set S, denote by 
y^g^M^ their disjoint uniorH. The categorical product of a finite family {Ms}ses of manifolds is 
the manifold 

\~\ Ms := {x : S ^ \_\ Ms \ x{s) G M, for ah s € S}. 

s£S ses 

We note that for each index s € 5" we have projection maps vr^ : Hs'gS ^-^s' ~^ are defined 

by 

•Ksix) = x{s). 

These projections are surjective submersions. 

We denote x(s) G Mg by Xg and think of it as s^^ "coordinate" of an element x G riseS-^*' 
Equivalently we may think of elements of the categorical product flseS"^'^'* ^ unordered tuples 
{xs)s&s with Xs^ Ms- 

^ The disjoint union may be defined by Usgs '■~ Usgs(^s ^ {•'})• 
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2.1.7 Remark. It is not hard to show that the categorical product, as defined above, has the 
fohowing universal property: given a manifold N and a family of smooth maps {fs : N Ms}s^s 
there is a unique map f : N \~\ses -^•^ with 

TTg o f = fg for all s € S. 

In fact categorical products are usually defined by this universal property [70] . 

2.1.8 Remark. For a family {Ms}s^s of manifolds indexed by a finite set S every ordering 
{si, . . . ,Sn} of elements of S identifies the categorical product flses-^-^'^* ^ manifold) with the 
Cartesian product Mg-^ x • • • x Mg^ . 

We are now in position to state: 

2.1.9 Definition (total phase space of a network {G,V)). For a pair [G^V) consisting of a graph 
G and a phase space function V we define the total phase space of the network {G,V) to be the 
manifold 

¥G = ¥{G,V) := |~| V{a), 

aeGo 

the categorical product of manifolds attached to the nodes of the graph G by the phase space 
function V. 

2.1.10 Example. Consider the graph 








a 



Define V : Gq ^ Man by ■P(a) = S*^ (the two sphere) and 'P{b) = S"^. Then the total phase space 
¥{G,V) is the Cartesian product 5^ x S^. 

2.1.11 Notation. If G = {0 ^ W}} is a graph with one node a and no arrows, we write G = {a}. 
Then for any phase space function V : Gq = {a} — > Man we abbreviate P({0 ^ {a\},V : {a} — > 
Man) = P({a},P : {a} -f Man) as Pa. 

2.1.12 Proposition. A map of networks ip : {G^V) — )• [G' ^V') naturally defines a map of corre- 
sponding total phase spaces 

¥ip : ¥G' -> ¥G. 

Proof. We use the universal property of the product PG = \~\a£Go^(^)'- define a map / from 
any manifold N to PG it is enough to define a family of maps {fa '■ N — > 'P{a)}aeGo- ^'^^ node 
a' of G' we have the canonical projection 

: PG' ^ V'ia). 

We therefore define 

(P(^), := : PG' ^ V'ivia)) = V{a) 
for all a € Go . □ 



9 



2.1.13 Example. Suppose G is a graph with two nodes a, b and no edges, G' is a graph with one 
node {c} and no edges, V'{c) is a manifold M, ip : G G' is the only possible map of graphs (it 
sends both nodes to c), and V : Gq ^ Man is given by V{a) = M = V{b) (so that V o (p = V. 
Then ¥G' ~ M, 



and P99 : M —^MxAI is the unique map with (¥'ip{x))a = x and {¥ip{x))b = x for all x G PG'. 
Thus : M — )■ M X M is the diagonal map x 1— > (rr, rr). 

2.1.14 Example. Let (G,P), (G',P') be as in Example EXE] above and -0 : {G' ,V') {G,V) 
be the map that sends the node c to a. Then P'i/' : PG — )■ PC is the map that sends (xq, Xb) to Xq. 

2.1.15 Remark. It is not hard to show that if a map of networks ip : {G, V) — )• (C, V') is surjective 
on vertices then ¥(p : ¥G' — >• PG is an embedding. If, on the other hand, (p is injective on vertices, 
then Fip is a surjective submersion. 

2.1.16 Remark. The total phase space map P : FinGraph/Man — t- Man is a contravariant functor: 
given two maps of networks 



To indicate that P reverses the direction of maps we now and subsequently write 

P : (FinGraph/Man)"" Man. 

The superscript °'' stands for the opposite category, i.e., the category with the same objects but all 
the arrows reversed. 

2.2 Open systems and their interconnections 

Having set up a consistent way of assigning phase spaces to graphs (that is, having set up networks 
of manifolds) we now take up a construction of continuous time dynamical systems compatible 
with the structure of the network. We build vector fields on total phase spaces of networks by 
interconnecting appropriate open systems. Our notion of interconnection is borrowed, to some 
extent, from the control theory literature. See, for example, Willems |71j . We therefore start by 
recalling a definition of an open (control) systems, which is essentially due to Brockett |72j . 

2.2.1 Definition. A continuous time control system (or an open system) on a manifold M is a 
surjective submersion p : Q ^ M and a smooth map F: Q ^ TM so that 



PG = {{Xa,Xi,) I Xa G V{a),Xb G V{b)} c^^MxM 



{G,V)^{G',V')^{G",.V") 



we have 



P(^ oip) = Fipo Ftp. 



(2.1.1) 



F{q) G rp(,)M 



Q 



F 



TM 



for all g G Q. (cf., for example, [73]). That is, the following diagram 




commutes 



M 



where vr : TM — > M is the canonical projection. 
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2.2.2. Given a manifold U of "control variables" we may consider control systems of the form 

F : M xU ^ TM. (2.2.1) 
Here the submersion p : M x U — > M is given by 

p{x, u) = X. 

The collection of all such control systems forms a vector space that we denote by Ctrl(M x [/ — > M): 
Ctrl(M xU ^ M) := {F : M xU ^ TM \ F{x, u) E T^M}. 

2.2.3 Notation (Space of sections of a vector bundle). Given a vector bundle E M we denote 
the space of sections of E ^ M hy TE or by T{E). 

Now suppose we are given a finite family {Fi : Mi xUi ^ TMi}fLi of control systems and we 
want to somehow interconnect them to obtain a closed system J^(-Fi, . . . , -F/v). This closed system 
is a vector field on the product Hi -^i- That is, J^{Fi, . . . , F/v) is a section of the tangent bundle 
r(P|- Mi)) \~\-Mi. What additional data do we need to define the interconnection map 

y : |~| Ctrl(Mi xUi^Mi)^r (r(|~| Mi)) ? 

i i 

The answer is given by the following proposition: 

2.2.4 Proposition. Given a family {pj : Mj x Uj — > Mj}jLi of projections and a family of smooth 
maps {sj : \~\Mi — )• Mj x Uj} so that the diagrams 



Pj 



commute for each index j , there is an interconnection map making the diagrams 

n, ctri(M, X f/, ^ M,) ^ r(r(n, Mi)) 

Ctrl(M, X U, ^ Mj) Ctrl(n, M, ^ M,) 

commute for each j. The components J^j of the interconnection map are defined by J'j[Fj) : = 
Fj o Sj for all j . 

Proof. The space of vector fields r(T(P|,. Mi)) on the product \~\ - Mi is the product of vector spaces 

r(r(nM,)) = rictri(riM, ^ m,). 

i j i 

In other words a vector field X on the product Hj is a tuple X = (Xi, . . . , X^), where 

Xj := D{prj) o X. 
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{D(prj) : T\~\Mi TMj) denotes the differential of the canonical projection pvj : \~\Mi — > Mj.) 
Each component Xj : \~\ - Mi TMi is a control system. 

To define a map from a vector space into a product of vector spaces it is enough to define a 
map into each of the factors. We have canonical projections 

TTj : |~~| Ctrl(Mi xUi^ Mi) Ctr\{Mj x Uj Mj), j = 1,...,N. 

i 

Consequently to define the interconnection map it is enough to define the maps 

: Ctrl(Mj X Uj ^ Mj) ^ Ctrl(|~~|Mi ^ Mj). 

for each index j. We therefore define the maps ^j : Ctr\{Mj X Uj Mj) Ctr\{\~\iMi ^ Mj), 
1 < j < iV, by 

y,{F,):=F,os,. 

□ 

2.2.5 Remark. It will be useful for us to remember that the canonical projections 

Wj : r(r|~|Mi)) ^ Ctrl(|~|M, ^ Mj) 

are given by 

Wj{X) = D{prj)oX, 

where D{prj) : T\~\ Mi — t- TMj are the differential of the canonical projections prj : \~\ Mi — )• Mj. 
2.3 Interconnections and graphs 

We next explain how networks of manifolds give rise to interconnection maps. To do this precisely 
it is useful to have a notion of input trees of a directed graph. Given a graph, an input tree I{a) 
of a vertex a is roughly, the vertex itself and all of the arrows leading into it. We want to think of 
this as a graph in its own right, as follows. 

2.3.1 Definition (Input tree). The input tree I{a) at a vertex a of a graph G is the graph with 
the set of vertices /(a)o given by 

/(a)o := {a}Ut-^{a)- 

where, as before, the set t~^{a) is the set of arrows in G with target a. The set of edges /(a.)i of 
the input tree is the set of pairs 

/(o)i := {(a, 7) I 7 e Gi, t(7) = a}, 

and the source and target maps I{a)i ^ /(a)o are defined by 

s(a,7)=7 and t(a, 7)=a. 

In pictures. 
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2.3.2 Example. Consider the graph G = (a 




(h^ as in Example 12.1.101 Then the 



a 



input tree I{a) is the graph with one node a and no edges: 



I{a) = 0. 



The input tree I{h) has three nodes and two edges: 




2.3.3 Remark. For each node a of a graph G we have a natural map of graphs 



i = ia-I{a)^G, ^„(a,7)=7. 



We stress that this need not be injective on nodes. For instance in Example 12.3.21 the map ^f, : 
I{b) G sends the nodes a and /3 to the same node a. 

2.3.4 Example. Consider the graph 



2.3.5 Remark. The input tree I[a) of a graph G is a directed tree of height 1: for any vertex x of 
/(a) with X ^ a there is exactly one edge with source x and target o. Also a is the only vertex of 
/(a) which is not a source of any edge (it's a root of 1(a)), and all the other vertices of I {a) cannot 
be targets of any edges (they are leaves of I (a)). 

2.3.6 Remark. It follows from Remark 12.3.51 above that if (p : I {a) — )• I{b) is an isomorphism of 
two input trees (these graphs may be input trees of two different graphs) then necessarily <p{a) = b 

2.3.7 Remark. Given a network (G, V) and a map of graphs (p : H G we get a map of networks 



a 




(2.3.1) 



The four input trees of G are the graphs 




ip:{H,Voy,) ^ {G,V) 



hence a map of manifolds 



Fif.FG^ ¥H. 



13 



2.3.8. Let {G, V) be a network and let a be a node of G. Consider the graph {a} with one node 
and no arrows. Denote the inclusion of {a} in G by La and the inclusion into its input tree I{a) by 
ja- Then the diagram 

{a}^I{a) 



G 



commutes. By Remarks 12.3.71 and 12.1.161 we have a commuting diagram of maps of manifolds 

P{a} P/(a) 



FG 



2.3.9. Let us now examine more closely the map ¥ja ■ P/(a) — )• Pa in 12.3.81 above. Since the set of 
nodes /(a)o of the input tree /(a) is the disjoint union 

/(a)o = {a}Ut^\a), 

and since ^0(7) is the source 5(7) for any 7 € t^^{a) C I{a)o, we have 

P/(a)=P(a)x |~| ^(5(7)). 

Since ja : {«} — ^ -^(a)o = {a} U t~^{a) is the inclusion, 

Pja : P/(a) ^ Pa 

is the projection 



Similarly 

is the projection 



V{a)x |~| V{s{-f))^V{a). 

7Gt-l(a) 

Fia : PG ^ Pa 



n ^ V{a). 

beGo 

Putting 12.3.81 and 12.3.91 together we get 
2.3.10 Proposition. For each node a of a network {G,V) the diagrams 

FI{a) = V{a) X fl^^t-Ha) ^(^7)) Via) 



commute. 
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2.3.11 Example. Suppose G = (a^ is a graph as in Example 1 2 . 1 . 1 1 and suppose 

a 

"P : Go Man is a phase space function. Then 

P/(6) =V{a) X V{a) x 

Pjb is the projection V{a) x V{a) x V{h) V{b) and 

Ctrl(P/(5) ^ P6) = Ctrl(P(a) x V{a) x ^ ^(6)). 

On the other hand P/(a) = V{a), Pja : 'P(a) 'P{(i-) is the identity map and 

Ctrl(P/(a) ^ Pa) = T{TV{a)), 

the space of sections of the tangent bundle TV {a), that is, the space of vector fields on the manifold 
V{a). 

2.4 Dependency 

Given a network (G, V) we have a product of vector spaces 

|~~| Ctrl(P/(a) ^ Pa). 

agGo 

The elements of the product are unordered tuples of {wa)a&Go '^^ control systems (cf. I2.1.6p . We 
may think of them as sections of the vector bundle 

Control(G,P) := |J Ctrl(P/(a) ^ Pa) ^ Go (2.4.1) 

over the vertices of G. This is the main reason for thinking of the collection of infinite dimensional 
vector spaces {Ctrl(P/(a) '¥a)}a&Go ^ vector bundle over a finite set. It will be convenient to 
have a notation for the space of sections of the bundle Control(G, P) — > Go. 

2.4.1 Definition. Let {G,V) be a network, as above. We refer to the bundle Control(G, "P) — )• Go 
as the control bundle of the network {G,V). We call the sections {wa)a£Go of the control bundle 
virtual vector fields on the network (G, P). We denote the space of sections by S{G,V). Thus 

g{G,V):= |~~| Ctrl(P/(a) ^Pa). 

aeGo 

We now argue that an application of the interconnection map cJ^ : S(G, "P) — )■ rr(P(G, P)) 
turns these "virtual vector fields" into actual vector fields on the total phase space P(G, P) of the 
network. Indeed observe that Propositions 12.2.41 and 12.3.101 give us 

2.4.2 Theorem. For a network (G, V) there exists a natural interconnection map 

J : |~~| Ctrl(P/(a) ^ Pa) ^ r(rPG)) 

aSGo 
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with 



for all nodes a G Gq. Here Wa '■ r(TPG) Ctrl(PGo — ^ Pa) are the projection maps defined by 
ma{X) = D{¥La) o X (q.v. Remark\MM- 



2.4.3 Example. Consider the graph G = (a 



b) as in Example 12.1.101 with a phase 



space function V : Gq Man. Then the vector field 

X = J^{wa,Wb) : V{a) x V{b) TV{a) x TV{b) 

is of the form 

X{x,y) = {wa{x),Wbix,x,y)) for ah G V{a) x V{b). 

If G = ( a ) Cb) j and P : Gq ^ Man is a phase space function, then 



{J{Wa,Wb, Wc)) {x, y, z) = {Wa{x),Wb{x, X, y), Wc{y, z)) 
for ah {wa,Wb,Wc) € S{G,V) and ah {x,y,z) e V{a) x x P(c)). 



3 Modularity 

3.1 Symmetry groupoid of a network 

In this section we give one possible version of what it means for some of the open subsystems in 
the tuple of the constituent subsystems {wa)aeGo G S{G,V) on a network {G,V) to be "the same." 
We start with a pair of examples. 

3.1.1 Example. Consider the graph 




from subsection II. 4i Choose a phase space function V : Gq = {1,2,3} — > Man with V{1) = 'P{2) = 
V{2>) = M for some manifold M. Then a typical tuple of open subsystems in S(G, V) that defines 
the dynamics on the network is a triple of the form (/i : M x M ^> TM, /2 : M x M ^ TM, /s : 
M X M ^ TM). It make sense to require that /i = /2 = /a- We can do it because the input trees 
/(I), /(2), 7(3) and the corresponding networks {I{i),V o ^j), 1 < i < 3, are all isomorphic (here, 
as before : I{i) — )■ G are the canonical maps, q.v. Remark |2. 3. 3p . 

3.1.2 Example. Consider the graph 
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Again define a phase function V by setting V{i) to be the same manifold M for all i. An element 
of S(G, V) is then of the form 

(/i : M TM, /2 : M ^ TM, h : M x M x M ^ TM, U : M x M x M ^ TM). 

Now it does not make sense to require that /s = /i but it does make sense to require that /i = /2 
and /a = /4 (!). Note that in this example the networks [I{1),V o and {I{2),V o ^2) are 
isomorphic as are the networks (1(3), 7^ o ^3) and (1(4), P o ^4). 

If we were to set V{1) = V{2) = ^(3) = M and V{A) = N ^ M, then an element of S(G, V) 
would be of the form 

(/i : M TM, /2 : M ^ TM, h : M x M x M ^ TM, f^-.MxMxN^ TN). 

In this case setting /i = /2 would make sense but setting /s = would not. And while (/(I), 'Po^^) 
and (1(2), V 0^2) would still be isomorphic, the networks (1(3),^ 0^3) and (/(4),Po^4) would not. 

3.1.3 Remark. In Example 13.1.11 there are 3^ isomorphisms 

ipi, : {I{j),Voij) ^ {I{i),Voii), 1 <i,j <3 

with 

for all i,j,k, and with 

(fa = id 

for all i (and consequently (pji = 'Pij^)- These 9 maps are an example of a groupoid. 
We recall the shortest definition of a groupoid: 

3.1.4 Definition. A groupoid is a category with every morphism an isomorphism. 

3.1.5 Remark. One may think of a groupoid IHI as a directed graph {Hi ^ Hq} together with an 
associative multiplication of pairs of edges with matched source and target: 

[a i— b i- c) [a i — b) 

an inversion map 

{a b) (a - — > b) 

and a unit edge ida : a a for every vertex a of H. We refer to the elements of Hq as the objects 
of the groupoid IHI and to the elements of Hi as isomorphisms of H. 

3.1.6 Example. In Remark 13.1.31 the groupoid G associated to the network {G,V) has three 
objects, namely the networks [I{i),V o ^j), 1 < i < 3, and 9 isomorphisms ipij, 1 < i, j < 3. For 
the corresponding graph see ()3.1.ip below. 

3.1.7 Definition (Symmetry groupoid G{G,V) of a network). The symmetry groupoid G = 
G(G, V) of a network (G, V) is a category with the following sets of objects and isomorphisms, 
respectively. The set of objects Go of G is the set of input networks 

{(/(a),Poea)WGo- 

The set of isomorphism Gi of G is the set of all possible isomorphisms of the input networks. 
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3.1.8 Example. Consider the network of Example I3.1.1[ As we have aheady pointed out the 
symmetry groupoid G of this network has 3 objects and 9 isomorphisms. It can be picture as 




(3.1.1) 



3.2 Groupoid-invariant vector fields 

Given a network (G, V) with a groupoid symmetry we should be able to talk about invariant 
elements of the vector space 'S>{G,V), the vector space of constituent open subsystems. This is 
indeed the case. There are several ways of making sense of invariants. The most concrete cuts 
out the subspace of invariant by appropriate equations. To set this up we need a number of short 
technical lemmas. We formulate them in a generality that is not needed immediately but will be 
useful later. The point of the lemmas is to prove that for a given a network (G, V) there is a natural 
action of its symmetry groupoid G{G,V) on the vector bundle Control(G, "P) — )■ Gq. 
We start by spelling out what we mean by the action of G{G,'P) on Control (G, "P). 

3.2.1 Notation. Denote the category of real vector spaces and linear maps by Vect. 

3.2.2 Definition. An action of the groupoid G{G,V) on the vector bundle Control(G,P) is a 

functor 

p : G{G,V) Vect, 

so that 

p{I{a),V o = Ctrl(P/(a) ^ Fa) 

for all nodes a of the graph G. Here as above Vect denotes the category of vector spaces and linear 
maps. 

3.2.3 Remark. The definition amounts to the following: 

1. For any two vertices a, 6 € Gq and an isomorphism ip : {I{a),V o ^q) — > (7(6), "P o ^h) in the 
groupoid G, there is an isomorphism 

p{ip) : Ctrl(P/(a) ^ Fa) Ctrl(P/(6) ^ P6); 

2. If {1(a), V o ^a) {1(b), V o ^fy) ^ {I{c),V o ^c) is a pair of isomorphism in G then 

°^) = p{'4') o p('/'); 

3. li ip : I{a) — )> I{a) is the identity isomorphism then p{<p) is the identity linear map. 

The notion of an action of a groupoid on a vector bundle is fairly old. See, for example, |74j . 
In the case where the vector bundle in question is a collection of vector spaces parameterized by a 
finite set of objects of the groupoid, as is the case for Control (G, P), it reduces to Definition 13.2.21 
above. 
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3.2.4 Lemma. Suppose ip : {G,V) — )• {G' ,V') is an isomorphism of networks. Then Ftp : PG' — > 
¥G is a diffeomorphism. 

Proof. Let ip"^ denote the inverse of iJj. Then ip o = id(G' -p) and o tp = id(G'/^-p/). Hence 

idpG = Fid(G,p) = r{ip o ^-1) = F{ip-^) o F^p. 

By the same argument 

idpG' =FipoF{'ip-^). 

Hence Fxp is invertible with the inverse given by F{ip~^). □ 

3.2.5 Remark. Here is a one hne category-theoretic proof of Lemma 13.2.41 since P is a functor, 
it takes isomorphisms to isomorphisms. 

3.2.6 Lemma. Suppose {G,V), {G' ,V') are two networks, '■ I{o) — > G the input tree of a vertex 
a of G, S^a' ■ G' the input tree of a vertex a' of G' and ip : I{a) — ?• I{a') an isomorphism of 
trees with "P' o o = "P o ^„ . Then the linear map 

Ctrl((^) : Ctrl(P/(a) ^ Pa) ^ Ctrl(P/(a') ^ Pa') (3.2.1) 

defined by 

Ctr\{ip) : (F : P/(a) TFa) ^ D{Fip\[^y)-'^ oFoFtp (3.2.2) 

is an isomorphism. Here f\{a} ■ {'^} ~^ restriction of ip to the subgraph {a} of G (by 

Remark \2. 3.6\ ip has to send a to a'); it is an isomorphism. 

Proof. By assumption (p : {I{a),'P o £,) ^ {I{a'),V' o is an isomorphism of networks. By 
Lemma [3.2.4l the maps Ff and Pi^lj^j are diffeomorphisms. Therefore Ctrl((^) has an inverse given 
by 

(F' : P/(a') ^ TPa') ^ D{F{ip\{a})) oFo Fip"'^. 

□ 

It fohows that we may define the functor p : G{G, V) — > Vect on isomorphisms of the groupoid 
G(G, V) by setting 

p{^) ■■= Ctrl(v;). 

3.2.7 Lemma. Given three networks {G,'P),{G' ,V') and {G",V"), and a pair of isomorphism of 
input networks 

iI{a),V o ^ {Iia'),V' o e)) ^ iIia"),V" o e") 

we have 

Ctrl(V' o 99) = CtrlV' o Ctrly?. (3.2.3) 
Proof For F e Ctrl(P/(a) Fa) we have 
Ctrliip o ip)F = F)(P(^o(^)|-J-p oFoP(^/>o(^) 

= Z)((P(^||(jj o P^||^/|)~"^) o F o Fip o Fip (since P is a contravariant functor) 
= D{Fi;\{,,})-') o (Z)(P(^|{„,})-1) o F o P(^) o P^ 
= Ctrl(^)(Ctrl(99)F). 

□ 
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We are now in the position to prove the main result of the section. 

3.2.8 Proposition. The symmetry groupoid G{G,V) of a network {G,V) acts on the vector bundle 
Control (G, "P) Gq. The action is given by 

[{I{a),V o ia) ^ {Hb),V o et)) ^ (^Ctrl(P/(a) ^ Pa) Ctrl(P/(6) ^ P6) 

where Ctrl (99) is defined by (j3.2.2p . 

Proof. We need to check that the three conditions hsted in Remark 13.2.31 hold for p{ip) = Ctrl (93). 
The first one holds by Lemma [3.2.61 The second by Lemma [3. 2. 71 Note finally that by construction 
if (p : 1(a) I (a) is the identity isomorphism then Ctrl (99) is the identity linear map. We conclude 
that the functor 

p= Ctrl : G{G,V) ^ Vect 
defines an action of the groupoid G{G,V) on the vector bundle Control (G, "P). □ 

Our next step is to define the space of invariant sections of the vector bundle Control(G, V) — > Go 
for this action, which is, by definition, the space of invariant virtual vector fields on the network. 

3.2.9 Definition (Invariant virtual vector fields on a network). Let (G, V) be a network We define 
the space V(G, V) of groupoid-invariant virtual vector fields on the network to be 

VG = V(G,P):= (3.2.4) 
{{wa) G S(G,P) j Ctr\{a)wa = for ah a £ G{G,V) with a: I{a) I{b)}. 

3.2.10 Example. Consider the network of Example 13.1.11 It is easy to see that 

V(G,P) = {(/l,/2,/3) G §(G,P) I /l = /2 = /s}, 

where, as before fi : M xM ^ TM are control systems. Note that in this case the space of invariant 
virtual vector fields is naturally isomorphic to the space Ctrl(P/(l) — >■ PI) = Ctrl(M x M — >■ TM). 
Note also that Ctrl(M x M ^ TM) is the space of invariant virtual vector field for the network 
{G',V) where 

G'= O 

is the graph with one vertex and one edge and the function V assigns the manifold M to the one 
vertex of G' . 

3.2.11 Remark. The reader may wonder in what sense the sections in V(G, P) are "invariant." 
There are several ways to answer this question. We start with the most concrete. Note that the 
space of H -invariant vectors for a representation p : H ^ GL{W) of a group H satisfies 

= {w£W\ p{a)w = w for all a G H}. (3.2.5) 

It is easy to see now that (|3.2.4p is an analogue of ()3.2.5p for groupoids. 

More abstractly, we note that the space is the limit of the functor p : Vect. Here H 
denotes the category with one object * and the set of morphisms Hom(*, *) = H. Similarly it is not 
hard to see that V(G, V) as defined above by equation ()3.2.4p together with the evident projections 
V(G,-P) Ctrl(P/(a) Pa) is the limit of the functor Ctrl : G{G,V) Vect. 
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3.2.12 Remark. We would like to think of the image ^{YG) of the space YG of invariant virtual 
vector fields under the interconnection map : '§i{G,'P) TT¥{G,V) as the space of "groupoid- 
invariant vector fields" on ¥{G,'P). Note that this is not literally correct since there is no natural 
action of the groupoid G{G, V) either on the the tangent bundle rP(G, V) or on the space of its 
sections. 

3.2.13 Remark. As we observed in Remark 12.3.51 the graph underlying the input tree network 
{1(a), V o of a network {G,V) is a directed tree of height 1. If : Ti ^ T2 is an isomorphism 
of trees of height 1, then (p necessarily sends the root rtTi of the first tree to the root rtT2 of the 
second tree. Hence if 93 : {Ti,Vi) — >• (72,7^2) is an isomorphism of networks and Ti,T2 are trees of 
height 1, it makes sense to define 

Ctr\{(p) : Ctrl(PTi PrtTi) Ctrl(Pr2 PrtTa) 

by a slight modification of ()3.2.2p : 

Ctr\{ip)F := DF ((^IrtTi)"^ o F oFip. (3.2.6) 

The proof of Proposition 13. 2. 8) is then easy to modify to show that that Ctrl is a well-defined functor 
from the groupoid of height 1 tree networks and their isomorphism to the category Vect of (not 
necessarily finite dimensional) real vector spaces and linear maps. 

3.3 An alternative notion of modularity 

Throughout the paper we take the point of view that a network is a directed graph G together with 
an assignment of a phase space to each vertex of G, that is, a pair 

{G,V : Go — )■ collection of phase spaces). 

Golubitsky, Stewart and their collaborators in their work on coupled cell networks additionally 
attach colors to edges of graphs. They require that maps of networks preserve the colors. In 
particular edges of input trees acquire colors from their canonical maps into the defining graphs, 
and symmetry group oids consist of color preserving isomorphisms. Thus from the point of view 
of Golubitsky et al. we work with monochromatic graphs. The results of this paper do have their 
colored analogues. The proofs, mutatis mutandis are the same. See [75] . 

4 Fibrations and invariant virtual vector fields 

We proved in Proposition 12.1.12] that a map of networks ip : {G,V) — t- (G'jV) defines a smooth 
map Fif : PG' — >■ FG between their total phase spaces (going in the opposite direction). The map 
ip, in general, does not induce a map between spaces of vector fields on the phase spaces PG and 
PG'. Nor does it induce a map between the spaces of virtual vector fields S{G,'P) and S{G','P'), 
let alone the spaces of group oid-invariant virtual vector fields YG and YG' . There is, however, a 
natural class of maps of networks that does. Following Boldi and Vigna [65] we call them fibrations. 
The notion of a graph fibrations is old. It arouse independently at different times in different areas 
of mathematics under different names. See [76] for a discussion. 

The goal of this section is to prove that a fibration of networks (p : {G,V) — t- [G' ,V') naturally 
defines a linear map ip* : Y{G',V') — s- Y{G,V). In the following sections we show that the maps 
ip* and Fip and the interconnection maps of the two networks are compatible in the best possible 
way. Consequently fibrations of networks give rise to maps of dynamical systems. 
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4.1 Fibrations 



4.1.1 Definition. A map (p : G —?■ G' oi directed graphs is a fibration if for any vertex a of G and 
any edge e' of G' ending at ip{a) there is a unique edge e of G ending at a with 93(e) = e'. 

A map of networks : {G,V) — s- (G'jV) is a fibration if the corresponding map of graphs 
if : G ^ G' is a fibration. 

4.1.2 Example. The map of graphs 

7' 



02) s 



5' 

sending the edge 7 to 7' and the edge 5 to 5' is a graph fibration. 

4.1.3 Example. All the maps of graphs in (|1.4.4p are graph fibrations. If we define the phase 
spaces functions on the three graphs by assigning to every node the same manifold M then the 
corresponding maps of networks are fibrations. 

4.1.4. Given any maps 99 : G — )> G' of graphs and a node a of G there is an induced map of input 
trees 

ipa : I{a) I{^{a)). 

On edges of I{a) the map is defined by 

V9(a,7) := ((^(a),¥5(7)) 
(cf . Definition 12.3. ip . Moreover the diagram of graphs 

I{a)—^I{p{a)) 



G ^G' 

commutes (the map '■ -^(a) ^ G from an input tree to the original graph is defined in Re- 
mark [2X3|). 

4.1.5 Lemma. If ip : G ^ G' is a graph fibration then the induced maps 
of input trees defined above are isomorphisms for all nodes a of G. 

Proof. Given an edge {ip{a),j') of I{ip{a)) there is a unique edge 7 of G with (p{'^) = 7' and ^(7) = a. 
Consequently 99^(0,7) = (v('^);7')- It follows that (fa is bijective on vertices and edges. □ 

4.1.6 Corollary. If a map of networks ip : (G, "P) — ?> {G\V') is a fibration then 

ipa : (/(a),Poea) ^ (/((^(a)),P'oe^(,)) 
is an isomorphism of networks. 

Proof. Follows immediately from Lemma 14.1.51 above and the definition of an isomorphism of net- 
works. □ 



22 



4.2 Maps between spaces of invariant virtual vector fields 

The goal of this subsection is to show that fibrations of networks send groupoid-invariant virtual 
vector fields to groupoid-invariant virtual vector fields. Namely we prove: 

4.2.1 Proposition. A fibration 

^:{G,V)^{G',V') 

of networks defines a linear map 

ip* : S{G',r') S{G,V) 

between spaces of sections of control bundles, that is, between spaces of virtual vector fields on the 
networks in question. 

Moreover if* maps the space N{G\V') of groupoid-invariant virtual vector fields to the space 
Y{G,V). 

Proof Recah (q.v. Definition EH]) that S(G,7') = flaeGo Ctrl(PJ(a) Pa). We define 
if* : |~~| Ctrl(P/(a ) ^ Pa') ^ |~| Ctrl(P/(a) ^ Pa) 

by 

for all a € Gq. Evidently <p* is linear. 

We now argue that invariant sections get mapped to invariant sections. Consider w' € Y{G', V'). 
Let a : {1(a), V o ^„) — (1(b), V o ^b) be an isomorphism in the groupoid G(G,V). Since is a 
fibration, the maps 

: {I(a),Vo^,) ^ (/((^(a)),P' oe^(,)) 

and 

^b--{I(b),VoCb)^(I(ip(b)),'P' °^^ib)) 

are isomorphisms. Therefore 

(fboaoip-^ : (I(ip(a)),V' o ^^(^)) (I(Lp(b)),V' o ^^(^j) 

is an isomorphism of networks, hence an isomorphism in the groupoid G(G',V'). Since w' is 
G{G',V') invariant by assumption, we have 

Ctrl(v36 oao (/J-^)w;^(„) = w'^^^^y 

Since Ctrl is a functor on networks of height 1 trees (cf. Remark 13.2. 13^ it respects compositions 
and takes inverses to inverses. Consequently 

Ctr\(ipb) o Ctrl(a) o Ctr\(ipa)-^w'^^^^ = w'^^^y 

Thus 

Qtr\(a)(ip*w')a = Ctrl((T) o Ct,\(ipa)-^w'^f^^^ = Ctv\(ipb)-^w'^^j^^ = (ip*w')b, 
which proves that (p*w' € Y(G,V). □ 
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4.2.2 Remark. The proof above shows that a fibration (p : {G^V) — >■ {G' ,V') also mduces a 
fully-faithful map of groupoids 

G(99) : G(G,P) ^G{G',V') 

which is given by 

4.2.3 Remark. Here is an alternative, more geometric, way to think of Proposition 14.2. 1] and its 
proof. The collection of maps 

{Ctr\{ipa) : Ctrl(P/(a) ^ Pa) ^ Ctrl(P/((^(a)) ^ P9j(a))} 

define a map of vector bundles 

^ : Control(G,P) ^ Control(G', 

which restricts to an isomorphism on each fiber. Hence Control(G, T') Gq is the pullback of 
Control (G', P') Gq. Consequently sections of Control (G', T'') — > Gq pull back to sections of 
Control (G, "P) — ?• Gq. Moreover the vector bundle map (f intertwines the actions of the groupoids 
G{G,V) and G{G'V' . Hence invariant sections pull back to invariant sections. 

4.2.4 Remark. In section [6] below we show that somewhat surprisingly the map 93* : Y{G','P') — )> 
V(G, V) of Proposition l4.2.1] is always surjective. We also characterize the kernel of ip* . In particular 
if (p : (G, V) — 7- (G', "P') is a quotient map (in the setting of coupled cell networks the fibers of such 
(p are equivalence classes of a balanced equivalence relation) then p* : Y{G','P') V(G,'P) is an 
isomorphism. 



4.3 Fibrations and maps of dynamical systems 



The goal of this subsection is to prove that fibrations of networks give rise to maps between 
dynamical systems. This is arguably the main result of the paper. Here is a precise statement: 

4.3.1 Theorem. Let p : (G, "P) (G',P') be a fibration of networks. Then for any groupoid- 
invariant virtual vector field w' S YG' the map ¥p : PG' PG intertwines the vector fields J^'{w') 
and ^{p*w'): 

D{¥p)oJ'{w') = J?{ip*w')o'¥ip. (4.3.1) 



Equivalently the diagram 



TPG' 

J"{w') 

PG'' 



■TFG 

"G 



(4.3.2) 



commutes. 

4.3.2 Remark. Note that by Proposition 14 . 2 . 1] since w' is groupoid invariant virtual vector field on 
the network (G',P'), the pullback p*w' is a groupoid-invariant virtual vector field on the network 
(G,P), i.e., p*w' eY{G,V). 



24 



Proof of Theorem \4.3.1 , Recall that the manifold PG is the product flaGGo Hence the tangent 
bundle bundle TFG is the product flaGGo ^I^'^- particular for each node a of the graph G, the 
canonical projection 

TFG T¥a 

is the differential of the map ■ FG — )• Pa. Here, as before, La ■ {a} ^ G is the canonical 
inclusion of graphs. By the universal property of products, two maps into TPG are equal if and 
only if all their components are equal. Therefore, in order to prove that (14.3.2P commutes it is 
enough to show that 

DFia o J^{ip*w') oFip = DFia o DFip o J^'{w') 
for all nodes a G Gq. By definition of the restriction of 99 : G — > G' to {a} ^ G, the diagram 

{a}^^{ip{a)} 



G 



G' 



(4.3.3) 



commutes. By the definition of the pullback map (p* and the interconnection maps , J^' the 
diagram 



TFa 




T¥ip{a) 



PG 




commutes as well. We now compute: 



DFia o j{ip*w') o : 



{^{ip*w'))a o Fip by definition of ^{ip*w')a 



DF {t^ia)0^\{a})o^'{w') 
DF{ipo La) o ^{w') 
DF{La)oDF(poJ^'{w'). 



(4.3.4) 



by K3J^ 
by definition of 
since F is a contravariant functor 
by 



^'iw')^(^a) 



□ 



4.3.3 Remark. In Lemma 15.1.11 b elow we show that surjective fibrations of networks give rise to 
embeddings of dynamical systems. Since balanced equivalence relations of the groupoid formalism 
of Golubitsky et al. |39l - l6l] define quotient networks, each balanced equivalence relation give rise 
to a surjective maps of graphs and hence to surjective fibration of networks in our sense. Thus a 
special case of Theorem 14.3.11 generalizes one direction of the groupoid formalism correspondence 
between invariant subspaces and balanced equivalence relations from ordinary differential equations 
to vector fields on manifolds. More specifically Theorem 14. 3. II is a generalization, to manifolds, of 
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Theorem 5.2 (direction (b)) of [61] and of Theorem 9.2 of [5^. We do not attempt to estabhsh the 
converse. More specificaUy, we do not attempt to characterize submanifolds of total phase spaces 
of networks that are preserved by ah groupoid invariant vector fields — we are only speaking to the 
"forward" direction. 

5 Dynamical consequence of Theorem 14.3.11 

In this section, we will discuss the implications of Theorem l4.3.1l Consider a fibration ip : (G, V) — ?> 
{G',V') of networks. Then 93 defines a map ^po : Gq ^ G'q from the set of vertices of the graph G 
to the set of vertices of the graph G' . In general ipo is neither injective nor surjective. However 
if a graph fibration (p : G ^ G' is surjective on vertices, it is automatically surjective on edges. 
Similarly if a graph fibration : G ^ G' is injective on vertices, then it is injective on edges as 
well. From now on we simply talk about injective and surjective graph fibrations. 

Next observe that a given fibration ip : (G, "P) — )■ {G\V') can always be factored as a map onto 
its image followed by the inclusion of the image: 



Hence any fibration can be factored as a surjection followed by an injection. We next analyze 
surjective and injective fibrations of networks. 

5.1 Surjective fibrations 

5.1.1 Lemma. Suppose ip : {G,V) {G' ,V') is a surjective fibration. Then P93 : FG' — ?> FG is an 

embedding whose image is a "poly diagonal" 



Proof. Assume first for simplicity that G' has only one vertex * and = M. Then for any 

vertex a of G we have 



PG' = M and PG = M*^" , where as before Go is the set of vertices of the graph G. By Proposi- 
tion [2Tl2] the map ¥ip : M ^ M'^o of the form 



{G,r)^{ip{G),V')^{G',V') 



A, 



{x e PG I 

— whenever (p{a) — (p(b)}. 



V{a) = V'{ip{a)) = V'{*) = M, 



¥ip{x) 



(x, . . . 



x) 



for all x G M. 
In general 



is the product of maps of the form 




X I— )> (x, . . . 



x). 



□ 
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5.1.2 Example. We consider the following surjective fibration (p : {G,V) {G',V') of networks. 
We take G to be the graph 






2n-2 












2n 





(5.1.1) 



with 2n vertices {n > 2). We choose a phase space function V that assigns a manifold M to all 
odd numbered vertices and a (different) manifold N to all even numbered vertices. We take G' to 
be the graph 

(5.1.2) 



ip{n) 




with two vertices and two arrows. We set V'{a) = M and 'P'{b) = N. We define the surjective 
fibration ip: G G' hy setting 

J a, n odd, 
1 6, n even 

The corresponding total phase space map ¥(p : M x N (M x N)"^ is given by the formula 

F{x,y) = {x,y,x,y, ...,x,y). 

The groupoid G{G',V') is trivial. Consequently Y{G',T") consists of a pair of control systems 
w'^ : M X N ^ TM and w'j^ : N x M ^ TN. They interconnect to define a vector field J^'{w') : 
M X N ^ TM X TN with 

J^'{w'){x,y) = {w'^{x,y),w'f,{y,x)) 

for all {x,y) e M X N = ¥G'. 

The groupoid G(G, V) is not trivial: all input networks corresponding to odd numbered vertices 
are uniquely isomorphic. That is, given the vertices 2A; + 1 and 2£ + 1, k ^ i, there is exactly one 
isomorphism : I{2k + 1) — ;> I{2£ + 1) as well il^ek = '4^ki~'^ '■ -^(2^ + 1) — ;> I{2k + 1). Analogous 
statement holds for input networks corresponding to even numbered vertices. 

The pullback map 

if* :Y{G',V') ^Y{G,V) 

is easily seen to be given by 

and y{ip*w') G rT{M x iV)" is given by 

J'{ip*w'){xi,yi,.. .,Xn,yn) = {w'a{xi,yi),wl{yi,Xi),.. . ,w'a{Xn,yn),w[{yn,Xn)). 

It is clear that ¥'ip{M x N) is an invariant submanifold of the vector field J'{ip*w')^ as should be 
expected in light of Theorem 14.3.11 

5.2 Injective fibrations 

Consider an injective fibration 99 : {G,V) — > {G',V') of networks. Lemma 15.2.11 below shows that 
the map ¥ip : ¥G' — )■ ¥G of total phase spaces is a surjective submersion. Combining this with 
Theorem 14 . 3 . II we see that for any group oid-invariant virtual vector field w' € Y{G',V') the map 

Fp : {FG',J^'{w')) {FG,J^{^*w')) 
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is a projection of dynamical systems. In particular for any singular point x of the vector field 
^{ip*w'), i.e., the point where the vector field is zero, the fiber Fip^^{x) is an invariant submanifold 
of the vector field J^{w'). 

Note also that since the map of graphs ip : G ^ G' \s injective, ip : G ^ '■p{G) is an isomorphism. 
Since 99 is also a graph fibration, there are no edges of G' with the source in G' \ 'p{G) and target in 
the image piG). Thus the image of 99 is a subsystem of G' that drives the dynamical system on G' . 
In other words the notion of an injective fibration makes precise the intuitive idea of a subsystem 
driving a larger network. 

5.2.1 Lemma. Suppose ip : (G, "P) {G' ,V') is an injective fibration. Then ¥ip : ¥G' FG is a 
surjective submersion. 

Proof. Since ip : G ^ G' is injective, the set of nodes G'q of G' can be partitioned as the disjoint 
union of the image p{Go), which is a copy of Go, and the complement. Hence 

FG' ~ |~| V{p{a)) X |~| V'{a') ~ PG x \~\ V'{a'). 

With respect to this identification of FG' with PG x Ua'^^piGo) '^'W) the map Fp : FG' FG is 
the projection 

PG x PI V'{a') FG. 

a'^^iGo) 

which is a surjective submersion. □ 

5.2.2 Example. Consider the injective graph fibration 

G' 




(5.2.1) 



Choose phase space functions VjV so that i : (G, "P) — )■ {G',V') is a map of networks. By the 
discussion above, for any choice of a groupoid invariant virtual vector field w' € Y{G','P') the 
dynamics in the subsystem (PG, J^{i*w')) drives the entire system (PG', J^{w')). This is intuitively 
clear from the graph (j5.2.ip since there are no "feedbacks" from vertices 4, . . . , 10 back into 1, 2, 3. 

5.3 General maps 

As we observed in the beginning of the section any fibration p : (G, P) — t- {G',V') can be factored 
as a surjection onto its image 

p:{G,V)^{p{G),V') 

followed by the inclusion 

i:{p{G),V')^{G',V'). 
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It follows from the two subsections above that for any groupoid invariant virtual vector field lu' G 
Y{G','P') the map of dynamical systems 

¥ip : {¥G',J^'{w')) {¥G,J^{f*w')) 

factors as a projection of dynamical systems 

¥i : {FG',,y'{w')) ^ {Fip{G),,y'{i*w')) 

followed by the embedding 

{Fip{G),y'{i*w')) ^ {FG,y{ip*w')). 
5.3.1 Example. Consider the graph fibration 



G G' 




from Example 14. 1.21 Choose a phase space function V' on G' and define V : Go Man by 
V{ai) = V{a2) = V'{a), V{h) = V'{b). Then if : {G,V) {G',V') is a fibration of networks. It 
factors as ^ 




6 Spaces of invariant virtual vector fields 

The purpose of this section is to characterize further and more precisely the space of groupoid- 
invariant virtual vector fields on a network (G, V) and to understand better the pullback maps 
^* : Y{G',V') {G,V) induced by fibrations of networks ip : {G,V) {G',V'). 



6.1 The space Y{G,V) as a product of spaces of fixed vectors 

It will be useful to introduce a bit more notation. 

6.1.1 Notation. Given a network (G, 'P) we have an evident bijection between the set Go of 
vertices of the graph G and the set Go = {{1(a), V o (,a)}aeGo of objects of the groupoid G{G,V). 
It will be convenient to identify the two sets: 

G{G,V)o = Go. 

6.1.2 Definition (Automorphism group). For a vertex a of a graph G, hence for an object of the 
symmetry groupoid G{G, V) of a network (G, V) we set 

Aut(a) := {ip : {1(a), V o (1(a), V o | is an isomorphism of networks } 

Clearly Aut(a) is a group under composition. We call it the automorphism group of the vertex a. 
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6.1.3 Remark. Note that Aut(a) is the collection of isomorphisms of the groupoid G(G,V) with 
source and target a. 

By construction Aut(a) acts on the vector space Ctrl(P/(a) — >• Pa), the fiber of the bundle 
Control(G,P) Go- We denote the space of fixed vectors by Ctrl(PI(a) Pa)^"*^") 

6.1.4 Remark. In general given an object a of a groupoid H = {Hi ^ Ho} we have a group 
Aut(a) consisting of isomorphism of H with source and target a. 

6.1.5 Definition. Given a groupoid H we say that two objects a and 6 of H are isomorphic if there 
is an isomorphism 7 of H with source a and target b. 

6.1.6 Remark. It follows easily from the definition of a groupoid that being isomorphic is an 
equivalence relation on the objects. We denote the collection of isomorphism classes of objects of 
a groupoid IHI by Hq/Hi and denote the isomorphism class of an object a by [a]. 

6.1.7 Lemma. Let {G,V) be a network. The space Y{G,V) of groupoid invariant virtual vector 
fields is isomorphic (as a vector space) to the product 

y Ctrl(P/(a) ^Pa)^"*("). 

[a]GGo/Gi 

Here as in Remark \6.1.3\ Ctrl(P/(a) — > Pa)^"*'-'*^ is the space of vectors fixed by the action o/Aut(a). 

Proof. Suppose w € Y{G,'P) is an invariant section of Control (G, "P) Gq. Then for any node a 
of Go and any automorphism ^ G Aut(a) we have 

Ctr\{ip)wa = Wa- 

Hence Wa G Ctrl(P/(a) Pa)^^'*("). 

If a and b are two isomorphic objects in the groupoid G{G,V) by way of ip : {I{a),V o ^g) — )• 
{I{b),Voi,) then 

Wh = Q\.r\{'4))wh- 

It follows that if we pick representatives ai, • • • , oat € Gq of the equivalence classes in Gq/Gi then 
the restriction map 

TV 

V(G,P) ^ □Ctrl(P/(ai) ^Pai)^"*("'), w ^ {wa,, . . . ,Wa^) 

i=l 

is an isomorphism of vector spaces. This proves the lemma. □ 

6.1.8 Example. Consider the network [G,V) where G is the graph 




® ® ® 
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and V assigns the same manifold M to each vertex of G. Then the input trees of G are all of the 
form 




and the corresponding input networks are all isomorphic. Moreover they have trivial automorphism 
groups. Consequently 

N{G,V) ~ Ctrl(M X M ^ TM). 
This is quite small compared to the space of all vector fields on the total phase space PG ~ M^^ . 

6.2 Maps between spaces of invariant virtual vector fields 

The goal of this subsection is to understand the pullback map ip* : Y{G',V') — )• Y{G,T') between 
groupoid-invariant virtual vector fields induced by a fibration ip : {G,V) — )• {G' ,V) of networks. 
We will see that p* is always surjective. To describe its kernel we need the following concept. 

6.2.1 Definition. Let p : {G,V) — > {G',V) be a fibration of networks. The essential image 

essim(/7 C G'q of (p consists of all the vertices a' G Gq so that there is an isomorphism 

of input networks for some vertex a of G. 

We say that ip is essentially surjective if essimc^ = Gq. 

6.2.2 Example. The map 

G' 




of networks is not surjective. But it is essentially surjective if 'P'(i) = 'P'(i) for all 1 < i < j < 10, 
i.e., if we assign the same manifold to all vertices of the graphs. 

6.2.3 Theorem. Let p : {G,V) (G',P') be a fibration of networks. Then ip* : V(G',P') 
V(G, V) is surjective. The kernel of ip* is the space 

keip* = {w' € N{G',V') \w'^, =Q for all a' G essim(/9}, 

where essim ip is the essential image of ip defined above. 

Proof. For w' G Y{G','P') the pullback ip*w' is zero if and only if the component {ip*w')a = for 
all a G Gq. Since {ip*w')a = Ctr\{ipa)~^Wip(^a) (q-v. proof of Proposition and since Ctr\{ipa)~^ 

is an isomorphism we conclude that 

ip*w' = <^ '^'^{a) fo'^ a G Gq. 

Finally note that an invariant section w' G Y{G','P') vanishes on the image of ip if and only if it 
vanishes on the essential image of ip. □ 



31 



6.2.4 Corollary. If ip : [G,V) — ?> {G',V') is an essentially surjective fibration of networks then 
if* : V(G','P') —7- Y{G,'P) is an isomorphism. In particular (p* is an isomorphism if ip is surjective. 

6.2.5 Example. Consider the map i of networks in Example I6.2.2I Since i is injective and essen- 
tially surjective the map i* : Y{G', V') V(G, V) is an isomorphism. Compare with Example l6.1.81 
Clearly the map i is very far from being surjective. 

6.2.6 Remark. As we pointed out in Remark 14.3.31 in the groupoid formalism of Golubitsky et 
al. the quotient maps defined by balanced equivalence relations are surjective. Hence the spaces of 
groupoid invariant vector fields on a network and on its quotient by a balanced equivalence relation 
are always isomorphic. 
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